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[1] DERIVE im Mathematik- und Physikunterricht, Gunter Scheu 
F. Diimmler Verlag, Bonn, Diimmlerbuch 4592, 1994 

[2] Mathematik am PC, Einfiihrung in DERIVE, Bernhard Kutzler 
Soft Warehouse - Hagenberg, 1994 

[3] Mathematics on the PC, Bernhard Kutzler 
Soft Warehouse - Hagenberg, 1994 


First North American DERIVE User Group Meeting 

Sunday, November 20, 1994, 9:00- 13:00 
Room Europe 10, Walt Disney World Dolphin 
Lake Buena Vista,, Florida, USA 

Report by Bernhard Kutzler 

The first North American DERIVE User Group Meeting was organized in the frame of the 
International Conference on Teaching Collegiate Mathematics (ICTCM) Nov. 17-20, 1994. 
Despite the early Sunday morning time at the end of the conference and despite of the temp- 
tation of being in the middle of Disney World with its many fun parks, approximately 30 peo- 
ple attended the meeting, among them Bert Waits, one of the two founders of the ICTCM 
conference series, several delegates from Europe and one from Australia. 

I had the great honor of chairing the meeting and started by reading Josef's welcome note. It 
was only half as charming as if Josef would have spoken himself, but people told me that just 
by hearing his words they almost could see him there. What a marvellous compliment from 
those who know Josef personally. 

David Stoutemyer, one of the two fathers of DERIVE, was the first speaker. He gave some 
insight in DERIVE version 3. His lecture was followed by presentations of Jeanette Palmiter 
(USA), David Sjostrand (Sweden), Lisa Townsley-Kulich (USA), Terence Etchells (UK), 
Robert Mayes (USA) and Bernhard Kutzler (Austria). 

We had many nice discussions following the lectures, met with new DERIVE friends, and 
once again enjoyed sampling Hawaiian macadamia nuts (which were sponsored by Soft 
Warehouse, Inc.). It was a very nice start for more North American DERIVE activities. 


Bernhard, many thanks for your nice report. I am jealous, not because of Disney World, but because of 
the macadamia nuts. There is a challenge for us in Europe. Dear German members (and members from 
the countries in the neighbourhood. I'm glad to announce the 


2. Deutsches DERIVE User Group Treffen 
im Rahmen der MNU 1995, Niirnberg (Ostern 1995) 

Mittwoch, 12. 4. 1995 14.00 - 16.00 Uhr 

Gaste und Interessierte sind naturlich herzlich willkommen. Falls Sie den Wunsch haben, einen Kurz- 
vortrag zu halten, uns etwas vorzuzeigen, Wiinsche oder Anregungen darzulegen, ware ich fiir eine 
Vorinformation sehr dankbar. Ich hoffe sehr, mit Ihrer Hilfe ein kurzes Programm zusammenstellen zu 
kb nn en. Ich danke jetzt schon Herm Wolfgang Propper fur die Organisation des Treffens. Naheres 
erfahren Sie im nachsten DNL (Marz 1 995). 
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Liebe DUG-Mitglieder, 

Anfangs dachte ich, es wurde schwierig werden, 
far jeden DNL einen “Letter of the Editor 44 zu 
schreiben. Ganz im Gegenteil, ich habe immer 
Probleme, all das, was ich Ihnen so nebenbei noch 
mitteilen will in einer Spalte unterzubringen. In 
diesem DNL mochte ich Sie auf die bereits ange- 
kiindigte Utility von Sergey Biryukow hinweisen. 
Die zahlreichen Files finden Sie ebenso wie die 
beiden Demos im Unterverzeichnis <LABEL>. Es 
muB eine ungeheure Arbeit gewesen sein, die Zei- 
chensatze zu entwerfen. Sergey und sein Team bitte 
vor den Vorhang. 

Fur mich ist es besonders erfreulich feststellen zu 
konnen, daB der DNL vielfach nicht nur rasch 
durchgeblattert wird, um spezielle „Rosinen“ fiir 
seinen Bedarf zu finden, sondern, daB sich viele 
Leser ernsthaft mit den angebotenen Themen aus- 
einandersetzen. Ein Ergebnis finden Sie im User 
Forum, zwei weitere werden erwahnt. Sie sind so 
umfangreich, daB sie einen eigenen Artikel wert 
sind. 

Nun muB ich ein Versaumnis nachholen. Leider 
habe ich bisher vergessen G. Scheu's zweites 
DERIVE Buch vorzustellen. Das wird jetzt geme 
mit einer Bitte um Entschuldigung nachgeholt. Ich 
mochte besonders darauf hinweisen, daB hier zahl- 
reiche physikalische Anwendungen behandelt wer- 
den. Zufrieden, Gunter? 

Beachten Sie bitte die beiliegenden Informationen. 
Beim DERIVE Journal mochte ich nochmals daran 
erinnem, daB Sie als DUG-Mitglied einen betracht- 
lichen Preisvorteil genieBen. 

Nun habe ich eine gute und eine schlechte Nach- 
richt fiir Sie. Zuerst die schlechte: wir werden den 
DUG-Mitgliedsbeitrag erhohen, und zwar auf 
S 340.- fur Osterreich, DM 52.-, bzw. BP 22.- fiir 
Europa und US$ 40.- fiir Ubersee. Aber jetzt die 
(hoffentlich) gute: der Umfang des DNL wird ii- 
berproportional auf mindestens 40 Seiten erweitert. 
Fassen Sie das auch als Bitte um weitere Beitrage 
auf (Sie sind aber recht fleiBig! !). Kein Beitrag geht 
verloren, und wenn der eine oder andere langere 
Zeit nicht erscheint, dann hat dies ausschlieBlich 
technische Ursachen (Platzbedarf, Themenzugeho- 
rigkeit, Ubersetzung, ....). Fiir 1995 habe ich eine 
Ausgabe mit Schwerpunkt Geometrie vorgesehen, 
eine andere soli sich mit dem Geschehen in oster- 
reichischen Schulklassen beschaftigen. 

Beachten Sie bitte die beiliegende Rechnung fiir 
den Mitgliedsbeitrag 1995. Sie erhalten automa- 
tisch ein Quittung mit dem DNL# 18. 

Ich wiinsche Ihnen ein frohes Fest und ein erfolg- 
reiches Neues Jahr 1995. 

Mit den besten GriiBen bis zum nachsten Mai 


Dear DUG Members, 

Producing the first issues of the DNL I thought that 
it would be difficult or impossible to write a “Letter 
of the Editor” for each DNL. On the contrary, Eve 
problems to bring all the additional facts I want to 
tell you into one column. In this DNL Ed like to 
draw your attention to Sergey Biryukov's utility. 
You can find the numerous files together with two 
demos in the folder <LABEL>. It must have been 
an enormous work to design and edit the font sets. 
Sergey and his team, in front of the curtain! 

It is a great pleasure for me to see that many of you 
do not only browse the DNL to find special “rai- 
sins” for their own purpose but are dealing seri- 
ously with the offered items. You can find one 
result in the User Forum, two others are so volumi- 
nous that they are worth an own contribution. 

I must admit that I have to make good something. I 
don't know why, but I have forgotten to place 
Gunter Scheu's 2 nd DERIVE book on our Book 
Shelf until now. I want to point especially to the 
fact that a couple of physics applications can be 
found in this book. Satisfied Gunter? 

Please take notice of the information enclosed. 
Concerning the DERIVE Journal I would like to 
remind you that as a DUG member you benefit 
from a considerable reduction. 

Now I have one good and one bad news for you. 
Let me take the bad one first: We will raise the 
DUG Membership fee up to AS 340 for Austria, 
DM 52 or BP 22 for Europe and US$ 40 for other 
destinations. But here is the - hopely - good news: 
the contents of the DNL will be extended overpro- 
portionally to at least 40 pages (instead of 30). 
Please take this fact as a call for further contribu- 
tions. (I must confirm that you have been very 
active and reliable until now.) 

No contribution - large or small - is getting lost. If 
one or another will not be published for a longer 
time then this is caused only by technical reasons 
(amount of space, subject not suited for the issue, 
translation, ...). For 1995 Eve in mind to dedicate 
one issue to geometry and another one shall deal 
mainly with activities in Austrian classrooms. 

Please pay attention to the enclosed invoice for the 
renewal of your membership. You will automati- 
cally receive a receipt together with DNL# 18. 

My wife Noor and I wish you all personally known 
or not a Merry Christmas and a successful New 
Year 1995. 

Sincerely yours 

( 
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The DER1VE-NE WSLE TIER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tents of 30 pages minimum. The goals of 
the D-N-L are to enable the exchange of 
experiences made with DERIVE as well as 
to create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone: 43-(0)2275/8207 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the 
D-N-L. It must be said, though, that non- 
English articles are very welcome nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
D-N-L. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE Newsletter will be. 


Preview: (Contributions for the next issues): 

Stability of systems of ODEs, Kozubik, SLO 
Prime Iterating Number Generators, Wild, UK 

Graphic Integration, Probability Theory, Linear Programming, Bohm, AUS 

DERIVE in Austrian Schools, some examples, Lechner, Voigt, Eisler a.o., AUS 

Tilgung fremderregter Schwingungen, Klingen, GER 

Continued Fractions and the Bessel Functions, Cordoba a.o., ESP 

Turtle Commands in DERIVE, Lechner, AUS 

DREIECK.MTH, Wadsack, AUS 

IMP Logo and Misguided Missiles, Sawada, HAWAII 

“Reverse" Discussion of Curves, Reichel, AUS 

Reichel - Klingen - Bohm - Splines, A thriathlon, AUS & GER 

3D Geometry, Reichel, AUS 

Parallel- and Central Projection, Bohm, AUS 

Conic Sections, Fuchs, AUS 

A Trick for Plotting 2D Plots, Roanes & Roanes, ESP 

Setif, France; Vermeylen, Belgium; Leinbach, USA 

Lymer, FRA; Baum, GER; Kayser, GER and others 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag.Josef Bohm 
Herstellung: Selbstverlag 
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Peter Baum, Kassel Germany 

Peter Baum dealt with Th. Weths “Lexicon of 
Curves 44 and sent an extensive completion together 
with a couple of TIF files of nice pictures. Mr 
Baum uses only the polar form and he yields 
interesting results.The paper is too large to be part 
of the User Forum, so it will be a contribution in 
1995. But I would like to include one of his plots: 



Message 3225: From SOFT WAREHOUSE to PUBLIC about INSTALLING DERIVE ON A LAN 

The compact DERIVE and DERIVE XM executable files are extremely small by today’s standards. Also neither 
program uses overlay files while mnning (i.e. once loaded, they make no further access to the hard or floppy 
disk). Both these facts make installing and mnning DERIVE and DERIVE XM on a Local Area Network 
(LAN) simple and straight-forward. 

The following describes how to install DERIVE on a LAN. The procedure for installing DERIVE XM is the 
same, except where noted: 

1 . Login with system supervisor privileges and create a directory named DERIVE on the LAN file server. If 
installing DERIVE XM, name the directory DERIVEXM to avoid a possible name conflict with an existing 
DERIVE directory. Give DERIVE users ’’read only” rights to this directory and give the system supervisor 
all rights (see your LAN manual for details). Note that it is NOT necessary to establish a command file se- 
arch path to this directory. 

2. Using the MS-DOS COPY or XCOPY commands, transfer all the files on the DERIVE or DERIVE XM 
distribution diskette to the new directory. 

3. DERIVE users often need to save mathematical expressions and/or graphics screen images in files. Nor- 
mally these files should NOT be stored in the DERIVE directory on the file server. Instead they should be 
saved in a working directory on a diskette, in a directory on the user’s own computer, or in the user’s own di- 
rectory on the file server. The DOS drive and directory current at the time DERIVE is started is the default 
working directory. Therefore, DERIVE should be started from the desired working directory and NOT from 
the DERIVE directory on the file server. 

4. On a LAN, it is usually preferable to start DERIVE using a DOS batch file. In addition to starting DERIVE, 
the batch file can automatically switch to the desired working directory and/or redirect printer output to the 
desired network printer. The batch file should be named DERIVE.BAT and saved on the file server in the 
standard batch file directory or any directory to which a command file search path has been established. 

For example, if the file server is drive F, the batch file 

0ECHO OFF 

F : \ DERIVE \ DERIVE . EXE %1 

starts DERIVE using the current DOS drive and directory as the working directory. If installing DERIVE 
XM, change the last line of the batch file to read 

F:\DERIVEXM\DERIVE.EXE %1 

Alternatively, the batch file 

@ECHO OFF 
C: 

CD C:\DERIVE 
F:\DERIVE\DERIVE.EXE %1 

makes the DERIVE directory on drive C the working directory. 

Finally, the batch file 
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@ECHO OFF 
F: 

CD F:\WORKING 
F : \ DERIVE \ DERIVE . EXE %1 

makes the WORKING directory on the file server the working directory. 


5. If you want to redirect DERIVE printer output to a network printer, you can include the appropriate printer 
redirection command before the DERIVE.EXE line in DERIVE.BAT (see your LAN manuals for details). 
For example, if you are running a Novell LAN and the Netware program CAPTURE.EXE is in the PUBLIC 
directory on drive F, the command line 

F:\PUBLIC\CAPTURE L=1 Q=name NT NFF TI=5 

redirects DERIVE printer output from the LPT1 port to the network printer named name. Also tabs are not 
expanded, extra form feed characters are not sent at the end of print jobs, and data is sent to the printer after 
A 5 second time-out. 

6. If the DERIVE user does not specify a drive and/or directory when attempting to load a utility file, DERIVE 
searches for the file in 

a. the DOS drive and directory current when DERIVE was started; 

b. the YDERIVE directory on the current drive; and finally 

c. the directory from which DERIVE was invoked. 

Thus, to provide access for all users to a utility file, copy the file to the DERIVE directory on the file server. 
Then the file will be found by search path c above. 

7. The directories listed above are also searched to find the DERIVE.INI initialization file when DERIVE is 
started, and to find the DERIVE.HLP help file when a request is made for help (see the DERIVE User Man- 
ual for details). 

8. You may want the state of DERIVE to be the same for all users when it starts. If this is other than the default 
initial state, save the desired state in a DERIVE.INI initialization file in the DERIVE directory on the file 
server (see Section 2.12 of the DERIVE User Manual for details). Then if individual users wish to modify 
the initial state, they can save their own DERIVE.INI files in their own working directory. 


Message 3226: From SOFT WAREHOUSE to PUBLIC about COMPUTING LARGE POWERS 

Normally DERIVE raises numbers to integer powers. For example, 2 A 3 simplifies to 8. However, before com- 
puting powers, DERIVE determines if the answer is going to be too big to store in memory. 

For DERIVE, the largest power it will compute is roughly 

2 A 50000 = 16 A 12500 - 3.16 10 A 15051. 

For DERIVE XM, the largest power it will compute is roughly 

2M00000 = 16 A 100000 ~ 9.96 10 A 120411. 

Although it does not raise numbers to huge powers directly, DERIVE can be used to convert such expressions 
to scientific notation using common logarithms like we learned our high trig class. 

For example, to convert 3 A 1000000, set the precision level to 12 digits, set the notation to decimal, and ap- 
proximate the expression 

LOG (3 A 1000000, 10) 

to give 477121.254719. Thus, 3 A 1000000 is approximately 10 A 0.254719 * 10 A 477121 
Finally, highlight and approximate the first factor to give 

1.7977 * 10 A 477121 which is 3 A 1000000 in scientific notation. 


Aloha, A1 Rich, SWH. 
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Message 3247: From SOFT WAREHOUSE to PUBLIC about 500,000 DIGITS OF PI 

In response to a challenge by a user, I had DERIVE XM compute n accurate to 500,000 digits on a 25mHz 486 
with 4Mb of memory. It took about a week to compute and another week to save the result in a file. 

It requires almost 50 pages to print using a very small type font. (No, I am not going to upload 500,000 digits on 
this BBS!) The last 20 of the 500,000 digits are 

4042487602513819524 

I would greatly appreciate someone independently (i.e. not using DERIVE XM) posting the last 40 of 500,000 
digits of 7i on this BBS so I can confirm that the DERIVE result is correct. Please include your source. 

Even on a 33mHz 486 with 8Mb of memory I was unable to get Maple V, Release 2 to compute more than 
about 100,000 digits of n before it exhausted memory. I would be interested to know if anyone can get Maple 
or Mathematica to compute 500,000 digits of n, and how much time and memory it takes. 

Aloha, 

A1 Rich, Soft Warehouse, Inc. 


Message 3250: From MICHAELWALSH to SOFT WAREHOUSE about #3247 / 500,000 DIGITS OF PI 

Hi Al— I've been involved in a discussion about computing tu on Prodigy. 

I know that tc/ 4 = 1 - 1/3 + 1/5 - 1/7 + 1/9 ... (this is a Taylor expansion of arcsin(l)) But this alternating se- 
ries doesn’t converge very rapidly. My question is: what algorithm did you use to compute 0.5M digits of tt? 
Thanks for your help. I’ve always been interested in number theory. Mike— Detroit 


Message 3255: From SOFT WAREHOUSE to MICHAELWALSH about #3250 / 500,000 DIGITS OF PI 

The algorithm DERIVE uses to compute n is based on Ramanujan's formula published in the Quarterly Journal 
of Pure and Applied Math, Volume 45, page 350, 1914. The formula is that 4/n equals the sum of the terms 

( -1 ) A m (1123 + 214 60*m) (1*3*5 ...* (2*m-l) ) (1*3*5* ...* (4*m-l) ) 


882 A (2*m+l) 32 A m (m!) A 3 

as m goes from 0 to infinity. DERIVE then uses some "tricks" to compute the series efficiently. 

I am still waiting for someone to post on this BBS the last 40 of 500,000 digits of re so I can verify DERIVE 
XM's result and sleep better at night. 

Aloha, 

Al Rich, Soft Warehouse, Inc. 

Message 3256:From SOFT WAREHOUSE to PUBLIC about MYSTERIOUS TRIG IDENTITY 

In response to a user inquiry, I stumbled across the following trig identity 

2 SIN ( 7i / 1 4 ) COS(7i/7) = COS (n/1) - 1/2 or equivalently 

SIN ( 3n/ 14 ) - SIN (7i/14) = COS(ti/7) - 1/2 

It looks like it should be easy to derive using the various rules for simplifying trig products, half-angles, etc. 
However, I have been stumped at finding a derivation. 

If there is a generalization of this identity, it could be incorporated into DERIVE's trig simplifier. Then expres- 
sions on the LHS of the identity could be transformed to the obviously simpler equivalent expressions on the 
RHS. 

However, the first step to finding such a generalization is to derive the above identity. I would greatly appreci- 
ate if someone could provide a derivation and/or generalization of this identity. 

Aloha, 

Al Rich, SWH. 
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Message 3263: From HADUD to SOFT WAREHOUSE about TRIG.IDENTITY 

Your formula can be derived from the trivial equality 

si "f< 7 - 3 >n) =c ° s (jf 


For brevity let E - e 14 . Then the LHS may be written 


E 4 _ E -4 _ (E + E~ l \E 3 - E + E~ l - E~ 3 ) 


2 i 


2 i 


cos — 

r-pi • 3 tc • tc 14 o 9 tc 3 2n 1 ^ 

Thus sin— -sin— = = 2cos 2 — = cos— — . Q.E.D. 

14 14 0 n 14 2 14 2 

2 cos — 

14 

I cannot believe that this is so clear for the average student. In my opinion it needs some knowledge, 
experience, skills and phantasy in manipulating trig expressions, Josef 

What is necessary to know? 


sin x cos y - sin(x + y) - sin(y - x) 


sm | y - x | = cos x 


( 1 ) 

( 2 ) 


pix _ p-ix pix p-ix 

sin x = and cos x = 

2 i 2 

(a 4 - b 4 ) = (a + b)(a 3 - a 2 b + ab 2 - b 3 ) 

cos(2x) = 2 cos 2 (x) - 1 

cos(3x) = 4 cos 3 x - 3 cos x 


(3) 

(4) 

(5) 

( 6 ) 


~ • n 2n • 3 7i • 7i 

2 sin — cos — = sin sin — 


14 14 


14 


14 


(i) 


. | 7t \ . (in 3 n\ . (n 3 n\ (3 n 

sm|<7-3)-j = s , n| i ---j = s ,„^---j = cos^- 


( 2 ) 


COS 




4 ni 4 ni 

gTT - £~YT _ E 4 - E~ 4 _(E + E~ l )(E 3 - E 2 • E~ l + E • E~ 2 - E ~ 3 ) 

2 i 2 i 2 i 


E 3 - E~ 3 E l - E~ 


2 i 


2 i 


C E + E - 1 ): 


f 3ni 

e 14 -e 

2 r 


3 ni 

IT 


e i4 — e 14 
2/ 


^i4 + e 14 



(3,4) 



3t7 a , 7T" 7T" 

cos-— 4 cos 3 — -3 cos— ~ 

14_ = 14 U = 2cos^-^ 

n ~ 7i 14 2 


2 ji 

COS — 

14 


2 ji 

COS — 

14 


1 


2n 1 


n 


1 


cos 2 — = cos — — - — = cos — 


14 


14 2 


7 2 


qued 


(5,6) 


We continue with Hadud's mail: 


More generally, consider 


sin — — KlIL - cos— with n, A: both odd. Letting a = — and proceeding as above we get 
2 n 2 n 2 n 


sin (( n-k - 1) a) - sin ((n — k — 3) a) H — ± sin a 


cos ( ha ) 
2 cos a 


plO 
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If k is fixed at 3 we have the identity 

sin ((ft - 4) a) - sin ((ft - 6) a) + sin ((ft - ±sin <2 = cos(2ft) - ^ with a - 

By interchanging the roles of sines and cosines corresponding formulas may be found for cosine terms on the 
LHS. 

Message 3265: From HARALD LANG to SWH about #3256 / MYSTERIOUS TRIG IDENTITY 

Perhaps the following simple identity can help: 


1/2 

6 -COS 


SIN 


14 


ATAN 


1 12 


r 

^ 7T ^ 

2 

" 7T ^ 

> 

- 108- COS 


+ 108- COS 


+ 5 


^ 7 ^ 


. 7 j 



1/2 



r 

" 7 T ^ 

> 


9- 

2- COS 


- 1 




. 7 > 

j 

j 


Message 3266: From HARALD LANG to SWH #3256 / MYSTERIOUS TRIG IDENTITY 

More seriously, simplifying trig expressions enters when DERIVE solves 3- and 4-degree equations. For exam- 
ple, try solving 

2 z 3 - z 2 - z + — = 0. 

4 

In my version of DERIVE, I get three roots, one of which is (the other two look similar) 


APPROX 


77 -COS 


ATAN 


73 


APPROX 


COS 

V v 7 jj 


= 0.9009688679 


0.9009688679 


But this is actually equal to COS(pi/7), a much simpler expression. This follows easily from the SWH/Hans 
Dudler formula. (I just figured out that I should use "Std Char set" and print to a File in derive; then I can 
download the formula just as is!) — Harald Lang 

DERIVE 6 produces another - but equivalent output for the solutions of the equation given above. 


SOLUTIONS 


77-SIN 


3 2 1 

2 ■ z - z - z + = 0 , z 

4 


f 

r yi 

ATAN 

— 


L 9 J 


3 ^ 


77 -SIN 


ATAN 


73 


3 6 6 

[0 , 2225209339 , -0 . 6234S9S01S, 0 . 90096SS679] 


77 -COS 


ACOT 


73 
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Message 3269: From JERRY GLYNN to PUBLIC about 500000 DIGITS OF PI 

A1 Rich recently calculated 500000 digits of pi on Derivexm. He asked for confirmation on this BBS. We are 
happy to announce his last 40 digits match up with two different sources so it look like congratulations are in 
order! ! ! 
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Message 3284: From HARRY SMITH to SOFT WAREHOUSE about 500,000 DIGITS OF PI 

Al, I know your 500000 digits of n has been confirmed, but I thought you might be interested in my program 
for computing ni. n computed accurate to 500,000 decimal digits has the following as its last 60 digits: 

969593993754953622322221974659619332529074042487602513819524 

The next 26 digits are: 2697391017563719753430045 

This was computed by the program: 

PiW - Compute PI to a million or so decimal places in 
Windows Version LOO, last revised: 1992/12/16, 0600 hours, 

Copyright (c) 1981-1992 by author: Harry J. Smith, 

19628 Via Monte Dr., Saratoga, CA 95070. 

All rights reserved. 

I used PiW to compute n to 500,000 decimal digits. This was done on an IBM AT compatible 33 MHz 486 DX 
computer using Windows 3.1 with 16 megabytes of RAM and 22 mega bytes of virtual memory. It took 37.3 
hours for algorithm (a) and 31.0 hours for algorithm (b). The results were the same. The divide and square root 
routines have been improved since these runs, so the program is a little faster now (10 to 20%). 

If I upload the program it will be in the two files PiWIOO.Zip and PiW500K.Zip. 

Its online description is: PiW vl .00 By Harry J. Smith, Computes Pi to a million decimal places 

Key words: Pi Math Windows C++ Precision FFT 

PiW is a program to compute n to a million or so decimal places on an IBM compatible PC using Microsoft 
Windows for memory management. With 16 megabytes of RAM it can compute n to 524,200 decimal places. 
With 32 megabytes of RAM it can compute n to 1,048,500 decimal places. 

PiW uses a Fast Fourier Transform (FFT) to speed up multiplications. It is written in Borland C++ for Windows. 
All source code is included. Also includes a copy of n to 500,000 places that was computed in 1.55 days using a 
33 MHz i486. It is a good example of object oriented programming (OOP). 

7i to 500,000 places is in file PiW500K.Zip. Every thing else is in file PiWIOO.Zip. — Harry 

Message 3288: From CHRIS LAMOUSIN to JERRY GLYNN about PLOTTING 

I'm having trouble plotting a function with Derive. I'm still figuring the program out. The function 

is y=x - (3/2)x A (2/3). I can get f(x) for x>0, but not for x<0. The program is set for all real numbers. Any Ideas? 

Message 3289: From JERRY GLYNN to CHRIS LAMOUSIN about #3288 / PLOTTING 

Yes I know what you need. Do Manage Branch Real and replot. Derive thinks correctly, that (-8) A (2/3) is a 
complex number. However this is an uncomfortable viewpoint in Calculus or before so Derive has this flexibil- 
ity. You'll probably want to leave the setting this way so do Transfer Save State. 

In commonly asked questions at the back of your manual you'll see this discussed. 

Message 3303: From JERRY GLYNN to PUBLIC about INTELLIGENCE? 

In all symbolic algebra programs I've tried 2 A 3001/2 A 3000 has produced a correct answer of 2 but has obviously 
done this by calculating both terms and dividing. We would not do it this way. All of these programs produce 2 
as an answer to 2 A (x+l)/2 A x so it's not that they can't do algebra but rather they can't do algebra when no vari- 
ables are present. Could this ability pass for an example of human intelligence? Are there other examples out 
there? If we knew what humans were good at and what algebra systems were good at maybe we could see how 
to work together and then maybe train people for these jobs. 

Do symbolic algebra systems simulate intelligence or are they just wonderfully fast calculators of symbolic 
expressions? 


Message 3304: From GREG SMITH to JERRY GLYNN about #3303 / INTELLIGENCE? 

The fact that you had to ask the question at all hints strongly at the answer ... 
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Message 3308: From HARALD LANG to JERRY GLYNN about INTELLIGENCE? 

I have a similar experience. It seems like DERIVE has a difficulty to handle redundant information. I have 
pointed out earlier the curious fact that she correctly simplifies 

1/3 1/3 1/3 

(x + 2) + (8-x + 16) - (27-x +54) =0 

to 0, but if you replace x by SQRT(5), then she fails! I don’t know if this is typical for computer algebra pro- 
grams in general, but no doubt we humans are not fooled by x being replaced by a number. — Harald 

Message 3309: From JERRY GLYNN to HARALD LANG about #3308 / INTELLIGENCE? 

Excellent example ... again algebra is something with variable(s) and arithmetic is something without variables 
and they are treated differently. 

When I wrote my comment I expected you to respond. How nice to be able to predict such a response at such a 
great distance. I’m inspired to think of more. 

Any more responses out there? 

Derive 6 works as expected: 

1/3 1/3 1/3 

0/5 + 2) + (8-^/5 + 16) - (27-^/B + 54) =0 

(In the next DNL you will find a fine example how an idea is growing on the Bulletin Board.) 


2D Plots Labeling 

Dr. Sergey V. Biryukov, Moscow 

Abstract 

Derive utility files for 2D plots labeling and axes drawing, numbering and labeling 
are described. Vector scaleable and rotatable fonts with -250 characters are avail- 
able. User fonts design is supported. 

Introduction 

Derive makes excellent 2D plots. It is for bright pupils, students and scientists. But ordinary pupils 
have some difficulties in keeping in mind plots and axes names, and in axes numbering calculation 
from scale values and cross position. Furthermore, axes cross is restricted to [0, 0] point. Our aim was 
to make 2D plots more clear and easy to percept. 

Utility overview 

Utility files set consists of: the main file label . mth, 5 font files f_* . mth, font designer 
label_fd.mth and appropriate document files *.doc. It supports plotting an arbitrary ASCII char- 
acter (except graphic and control ones) or sequence of characters (text) in the 2D Plot window. Char- 
acters position, width, height, rotation (slant) and characters step in a text can be defined. Axes with 
arbitrary point of their cross, numbers near ticks and labels near arrows can be plotted. 

One character plotting 

CHAR_(xy , wh , ch , rot) is the base function of the utility. It prepares a set of points for character ch 
plotting in Options>Display>Points>Connected and Small mode, xy = [. x_coordinate , y_coordinate ] - 
lower-left point of the character rectangle, and wh = [width, height ] defines rectangle width and height 
of this rectangle. Optional scalar parameter rot defines character rotation counterclockwise in radians. 
It is zero by default. It is necessary to preload a vector font (see next paragraph) to plot characters. 
CHAR_ ([1,2], [3, 4], "A" , 0) , for example, simplifies to 

[[4,2], [3. 5, 3. 5], [2. 5, 6], [1.5, 3. 5], [1,2], [1.5, 3. 5], [3. 5, 3. 5]] 

and gives a 3 units width and 4 units height upper case letter "A". The lower-left comer of the letter is 
at the point [1,2]. 
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prepares points for plotting num_x & numjy vectors of text under ticks. 

If the tick is in the point of axes cross appropriate text is shifted left or down. 
num_x= [text_for_left_tick, text_for_2nd_tick, . . .]. 

Three previous functions are combined in ALI AXES() one that has more clear input parameters and 

prepares points for plotting axes, axes labels and axes numbering. wh_x & wh_y in the second and 
third arguments are different and define characters width and height for the right neighbour text. 

ALL i AXES (axs, numbers, labels) = 


ALL AXES 


- 

xyl 



- 


xy2 

r wh x num x -i 

r Wh X 

label x -| 


axes_cross 

r \ \r 

L wh y num y J 

L wh y 

label_y -1 

- 

scale xy 



- 


Example 

#1 : L0 AD ( D : \ Df D\ DN L \ DN L94 \MTH16 \ L AB E L \ F_S , MTH ) 

#2 : LQAD(D: \DfD\DNL\DNL94\MTH16\LABEL\LABEL .MTH) 


#3: numb_x := [[-, 8], [- , 6], [-, 4], [- , 2], [0], [2], [4], [6]] 


#4: 

numb_y := [[- , 

6] 

, [-, 

4], [- 

, 2], [], 

[2], 

[4]] 







S' 


-S.5 

-6.2 ' 













S.B 

6.2 


1 1 

o 

UJ 

0.4] 

numb_x 


' [0.5, 

0.8] 

[X] ' 

#B : 

aa := ALL_AXE5 



0 

0 

1 

- [0.3, 

0.4] 

numb_y . 

1 

- [0.5, 

0.8] 

[Y] . 



‘v 


2 

2 . 










s' 













#6: 

la LABEL [1.7, 

1.2] 

, [0.4 , 

0.6], [0.5, 0] 

, [Y, =, 

x, 

o, 

., 7, 

5], 0, 



#7: [aa, la, y = x - 0.75] 
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(“DUG for YOU” and the arrows are added by me, then I tried to create a New Year's screen, Josef) 


#8 : L0AD(D : \ Df D\ DN L \ DN L94 \MTH16 \ L AB E L \ F_E . MTH ) 

#9: LABEL ([-8 , 3], [0.5, 1], [0.6, 0.1], [D, U, G, , f, o, r, , Y, 0, U]) 

#10: SHIFT_LABEL(LABEL([-&, 3], [0.5, 1], [0.6, 0.1], [D, U, G, , f, o, r, , Y, 0, U]), [0.1, 0.1]) 

#11: [ARR0W([0, 0], [-7, 2.8], 0.5), ARR0W([0, 0], [-4.7, 3.2], 0.5), ARR0W([0, 0], [-2.3, 3.7], 0.5)] 

#12 : LOADED: \DfD\DNL\DNL94\MTH16\LABEL\F_E .MTH) 

#13: LABEL ([-3 . 5 , 1.5], [0.6, 1], [0.8, 0], [A, , H, a, p, p, y]) 

#14: LABEL ([-3 , -0.5], [0.8, 1], [1, 0], [N, a, w, , Y, a, a, r]) 

#15 : LOADED: \DfD\DNL\DNL94\MTH16\LABEL\F_NUMB3 .MTH) 

#16: LABEL ([-1. 5 , -2.5], [1.2, 1], [1.5, 0], [2, 0, 0, 8]) 


#17: LABEL [-4, -2.5], [0.3, 0.4], [0.3, 0], [2, 0, 0, 8], — , — 

L 10 3 J 


#1B: 

LABEL^[2.5, 2], [0.2, 0.3] 

1 1 

□ 

LU 

□ 

[2, 

0, 0, 

8], 

7 r 

5 




#19 : 

SHIFTJ.ABEL (LABEL ([-1.5, - 

2.5], [1.2, 

1], 

[1.5, 

0], 

[2, 

0, 0, 8]), 

[0.05, 

0]) 

#20: 

SHIFTJ.ABEL (LABEL ([-1.5, - 

2.5], [1.2, 

1], 

[1.5, 

0], 

[2, 

0, 0, 8]), 

[-0.05 

, 0]) 


#21: B0X([-5 , -3.5], [5.5, 3.5]) 


As you can see in the following figures, Sergey's tools are working with DERIVE 6 in the 
same way. 
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See first the font sets from 1994 and then the changed sets adapted for DERIVE 6. 
f_mef . mth : [ 'V , >" , " £ " , "n" , "=" , "±" , ">" , "<" , " f" , "J " , V , 'V ] 
f _mf 2 . mth : ° , V" , " n " , " 2 " , ] 

f _mf r . mth : [ ">" , "<" , " f" , "J " , 'V , 'V , " ° ' , ' 7 " , " " " , " 2 " , "." ] 
f nO . mth : [ "#" , "$" , , " [ " , "\" , " ] " , " A " "I" ] 


#1 : LOADCD: \DfD\DNL\DNL94\MTH16\LABEL\LABEL . MTH) 

# 2 : LOADCD: \DfD\DNL\DNL94\MTH16\LABEL\F_MEF . MTH) 

#3: LABEL([-4, 0,5], [0,5, 0.S], [0.7, 0], [®, avg_, e, n, ~i d_ , +, >, <, int_u, nnt_l, cHv_, appr_]) 

#4: LABEL([-2, -1], [0,3, 0,5], [0.4, 0,05], [®, avg_, e, n, “i d_ , ±, >, <, nnt_u, int_], div_, appr_]) 


00 0 € n = ±: 

1 \ a ° 

-4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 

oo j/j (C n _ 

0.5 L 1. 5 1 2-3- 2fS 3 3.5 4 

4e.> < f J ° 

— i 


#5 : L0AD(D: \ Df D\ DN L \ DN L94 \MT H16 \ LAB E L \ F_M f 2 . MTH) 

#6: LABEL([-2, 0.5], [0.25, 0.4], [0.3, 0], [°, •, ■, 7, n_, 2_, box_]) 


0 o - 

-0.5 

-2 -1.5 -1 -0.5 

0.5 


# 7 : LOADCD: \DfD\DNL\DNL94\MTH16\LABEL\F_N0 . MTH) 

#&: LABEL C [-4, 0.5], [0.5, 0.S], [0.7, 0], [#, S, @, [, \, ], A , ', {, |, }, -]) 


[ \ ] 

lA ' { } ~ 

-4 -3 -2 -1 

12 3 4 

— 1 


The next sets are still valid. 

f s.mth: "double_quotes" , ", ", "/", 

n on ii -i ii ii o ii no ii ii A ii ii uz n n c n n 7 n n on iiq n n . n n . n n ^ n n _ n n ^n n 9 n n 

^r^r^r-^r^-r^r^r ' r 0 r J r • r r r ^ r r ' r r 

"A", "B", "C", "D", "E", "F", "G" , "H", "I", "J" , "K", "L", "M", "N", "O", "P", "Q", 
"R" , "S" , "T", "U", "V", "W", "X", "Y", "Z"] 

f _t . mth : [" ! ", "double_quotes" , " 1 ", " (", " ) ", ", ", " . ", 

n o n n a n n on no n n a n n c: n n r n n 7 n n on iiq ii ii .11 n . n n ^ n n _ n n ^n n 9 n 

I j O j J j . , , , \ ^ 1 ' 1 r 

n n n ^ n n y n 11 Z 11 ] 


It is important that the two lists (vectors) forming the font-set-mth-file are called char_names 
(the short one) and chars (the huge one). Have a look at this vector and admire the enor- 
mous work lying in them. Maybe you will create another font set for us!! Editor. 
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Der Fermat-Punkt im Dreieck 
The Fermat-Point in the Triangle 

ELR.Geyer, Wiesbaden, GER 

Wahrend einer Projektwoche an unserer Schule (Gymnasium) im Fruhjahr 1993 
beschaftigte sich eine Gruppe von 7 Schiilem der Jahrgange 9 bis 13 mit den 
Einsatzmoglichkeiten von DERIVE im Mathematikunterricht. Bis auf einen Schil- 
ler hatten sie keine Erfahrungen mit dem Programm. Nach dem Einfiihrungstag 
(4 Schulstunden) hatten sich alle soweit in die Bedienung eingearbeitet, dass sie in 
der Lage waren, sich an aktuellen Problemen des Mathematikunterichts der jewei- 
ligen Jahrgangsstufe zu versuchen. 

Am 3. Tag stellte ich alien gemeinsam eine Aufgabe: 

Es ist die Lage desjenigen Punktes in einem (spitzwinkeligen) Dreieck zu suchen, 
dessen Abstandssumme zu den drei Eckpunkten minimal ist - der Fermat - Punkt. 

During a project week at our school (grammar school) in spring 1993 a group of 7 students 
(age 15 - 19) dealt with DERIVE. After a short introduction - 4 hours - they were able to use 
DERIVE in connection with actual problems of their curriculum. On the 3rd day I set them all 
the same task: 

Find the point in an acute triangie ABC with a minimal sum of distances to the points 
ABC - the Fermat point. To explain the problem I showed the following proof by construction:) 

Zur Erklarung fiihrte ich den folgenden geometrischen Konstruktionsbeweis vor: 


Sei T zunachst ein beliebiger Punkt im 
Dreieck AABC. Dreht man das Teildrei- 
eck ACAT um 60° in positiver Drehrich- 
tung um den Punkt A, so ergibt sich ein 
Polygonzug CTTB, dessen Lange die 
gesuchte Abstandssumme ist. Die Lange 
des Polygonzuges ist minimal, wenn er 
auf der Geraden C'B liegt. 

Analog dreht man das Teildreieck ATBC 
um B um -60°. Der gesuchte Fermat- 
Punkt F ergibt sich also als Schnittpunkt 
der beiden Geraden C'B und C "A. 



This is the original screen shot from 1994. The labels were 
produced using Sergey Biryukov's labeling tool. 


T is an arbitrary point in AABC. Having rotated ACAT about 60° in positive direction you ob- 
tain the sum of the distances with the polygon C'T'TB. Ist length is minmal, if it coincides 
with the line C'B. In the same way you can rotate ATBC round B by -60° (negative). The 
Fermat - point is the intersection point of C'B with C" A. 

It was interesting to observe how the pupils found different ways to solve the problem ac- 
cording to their mathematics knowledge. 
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We can use dynamic geometry to illustrate the problem and to verify that interection point FP 
is indeed the solution point: 


I show the Cabri-figure (PC and TI-92/Voyage 200) and the T\-Nspire-f\gure(s), Josef. 



Zur analytischen Bestimmung von F wurde vereinbart, das Dreieck in einem kartesischen Koordina- 
tensystem mit A(0/0) zu beschreiben. 

Die Schuler kamen nun selbstandig je nach Kenntnisstand auf folgende Losungsansatze: 

I. Berechnung des numerischen Abstandswertes als Funktion des Punktes F(xj;) und 
Minimieren durch sukzessive Einengung der Deflnitionsbereiche fur x und y 

Calculation of the numerical value as a function of the point F(x,y) and minimization by a 
step by step restriction of the domains for x and y 

Definition of the TOTAL sum of distances: 

#1: [bl:=, b2:=, cl:=,c2:=] 

2 2 

#2: da(x,y) := V(x +y ) 

2 2 

#3: db(x,y,bl,b2) := ^/C Cx— bl) +(y - b2) ) 

2 2 

#4: dc(x , y , cl , c2) := V(Cx-cl) +(y-c2) ) 

#5: T0TAL(x,y,bl,b2,cl,c2) := da(x,y)+db(x,y,bl,b2)+dc(x,y,cl,c2) 
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Substitution of the co-ordinates of the given vertices of the triangle: 

#6: L(x , y) := TOTAL(x , y, 7, 2, 3, 5) 

It follows an evaluation of function L in a 5 x 5 matrix, 
with increment h starting with x = s and y = t: 

#7: FERM(s , t, h) := VECTOR(VECTOR(L(x , y), x, s, s + 4-h, h) , y, t, t 

+ 4-h, h) 

#8: FERM(3 .41, 3.12, 0.002) 

' 10.30679813 10.30679185 10.30678803 10.30678668 10.3067878 

10.30679055 10.30678457 10.30678106 10.30678002 1.30678146 

#9: 10.3067845 10.30677882 10.30677562 10.3067749 10.30677665 

10.30677998 10.30677461 10.30677173 10.30677131 10.30677338 

. 10.306777 10.30677194 10.30676936 10.30676927 10.30677165 

#10: FERM(3 . 415 , 3.127, 0.001) 

' 10.30676992 10.3067701 10.3067709 10.30677232 10.30677436 

10.30676901 10.30676927 10.30677015 10.30677165 10.30677377 

#11: 10.30676848 10.30676882 10.30676978 10.30677136 10.30677356 

10.30676834 10.30676876 10.3067698 10.30677146 10.30677373 

. 10.30676858 10.30676908 10.3067702 10.30677194 10.3067743 

Value of the distance-sum-function: 

#12: L (3 . 415 , 3.13) = 10.30676834 

Die beiden Beispiele aus dem DERIVE-Programm zeigen, dab die Auswertung der Matrizen zwar 
schwierig - weil unubersichtlich - ist, aber die Lage des Fermat-Punktes laBt sich so doch gut einen- 
gen. 

Both examples show that the evaluation of the matrices is difficult because not easy to sur- 
vey, but we can find the position of the Fermat point approximatively. 

The numerical approach would be easier to follow if we could include the x- and y-values as well in the 
table. This is possible with little manipulating the matrix - and the respective function could be pro- 
vided as a tool by the teacher. 

Producing a table including the x- and y-values: 

#13: val ues(s , t , h) := APPEND([APPEND([ ] , VECTOR (y , y , t , t+4h , h)) ] , FERM(s,t,h) 1 ) 1 

#14: val ues(3 . 41 , 3.12, 0.002) 

3.41 3.412 3.414 3.416 3.418 

3.12 10.30679813 10.30679185 10.30678803 10.30678668 10.3067878 

3.122 10.30679055 10.30678457 10.30678106 10.30678002 10.30678146 

#15: 

3.124 10.3067845 10.30677882 10.30677562 10.3067749 10.30678146 

3.126 10.30677998 10.30677461 10.30677173 10.30677131 10.30677338 

. 3.128 10.306777 10.30677194 10.30676936 10.30676927 10.30677165 


Then proceed as shown above. 
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II. Differentiation der Abstandssummenfunktion 

(Differentiation of the distance-sum-function) 

Fur die Schuler der Klasse 1 1 lag es nahe, den Extrempunkt der Abstandssummenfunktion L(x,y ) iiber 
die Ableitung zu suchen. Als Funktion zweier Veranderlicher entzieht sie sich allerdings dem fur 
Schuler bekannten Losungsverfahren. Einer der Schuler hatte hier die Idee, die partiellen Ableitungen 
von L(x,y) als Funktion der einen Variablen bei konstantem Wert der zweiten Variablen zu plotten. 
Hierzu defmierte er einen Vektor, dessen Komponenten jeweils die so bestimmten partiellen Ableitun- 
gen waren. Die dargestellten Graphen fiihren zwar nicht direkt zu den Koordinaten des gesuchten 
Punktes, zeigen aber „auffalliges“ Verhalten in der Nahe der Koordinaten des Fermat-Punktes und 
insbesondere an den Grenzen der durch die Dreieckspunkte bestimmten Definitionsbereiche. 

The pupils of form 1 1 tried to find the extremal pont of the distance sum function using the derivative. 
Being a function of two variables it exceeds the pupils' known algorithms. One of the students had the 
idea to plot the partial derivatives of L(x,y) wrt one variable keeping the other one constant. The 
graphs don't lead directly to the point we are looking for ; but they are showing a very “strange“ behav- 
iour in the neighbourhood of the point and especially at the borders of the domain defined by the an- 
gles of the given triangle. 


Partial derivative of L(x,y) wrtx, which will be evaluated for 2 < y < 4, 
with an increment of 0.5 


VECTOR 


— LO, y), y, 2, 4, 


0.5 

j 


Partial derivative of L(x,y) wrt y, which will be evaluated for 3 < x < 3.5, 
with an increment of 0.1 


{ d 


VECTOR 


.dy 


LO, y), x, 3, 


3.5, 0.1 



Original figures from 


Caution: The earlier Derive versions were 
not able to plot implicit functions, and they 
accepted all variables as the independent 
variable. Now we can do implicit plots and 
the graph of the 2 nd family has another ap- 
pearance: 



DNL#16 (DOS-Derive) 
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II. Nachbildung des geometrischen Beweises als Schnittproblem zweier Geraden 
(Reproducing the geometric proof as intersection point of two lines) 

Hierzu muBten die Koordinaten der gedrehten Punkte C' und C" bestimmt werden. Die Schuler expe- 
rimentierten dazu zunachst mit Tangens und Pythagoras. Zur Vereinfachung wurden dann die erfor- 
derlichen Drehungen nur auf Drehungen um den Ursprung (Punkt A) beschrankt. AuBerdem stellte ich 
den Schulern die Matrix einer Drehung um den Ursprung und ihre Anwendung auf Ortsvektoren - 
sozusagen als “Black Box-Funktion“ zur Verfugung. Entscheidend war, dass die Koordinaten der ge- 
drehten Punkte C_ und B_ zur Verfugung standen. 

Damit konnten die Geradengleichungen aufgestellt werden und das Schnittproblem auf dem Niveau 
der Klasse 9 gelost werden. 

The co-ordinates of the roteted points C'and C“ were necessary. At first the students experi- 
mented with Tangens and Pythagoras. As a simplification I made the rotation matrix avail- 
able to the students - as a Black Box. Having the points C'and B' the problem was easy to 
solve for pupils even of level form 9. 

A ABC given by A(0,0), B(bx.by), C(cx,cy) 

y2 - yl 

#1: "hne(xl, yl, x2, y2, x) := (x - xl) + yl 

x2 - xl 

"ax ay 
bx by 

#2: tn'ang"le(ax , ay, bx , by, cx , cy) 

cx cy 
.ax ay . 

Rotation matrix around the origin = A(0,0) 

' C0S(cO - SIN (a) ‘ 

#3: rot(oO := 

. SIN (a) CO 5 (a) . 

Calculation of rotated points C' and B' (c_ and b_): 

#4: c_(cx, cy) := notf l-[oc, cy] 



#6: b_(bx, by) rot| - J-[bx, by] 

I" bx ■ by by 4 3 ■ bx 

#7: b_(bx , by) := + , - 


Line 

s connec 

ting 

B with C 

(line G) ar 

id B 1 with C (lira 

3H): 





G(x , bx 

. by 

. cx, cy) 

:= "h'ne((c. 

— (cx , cy)) , (c. 

_(cx, cy 

)) , bx, by) 



#8: 





1 


2 




H(x , bx 

i by 

i cx, cy) 

:= "h'ne((b. 

-(bx, by)) , (b. 

_(bx, by 

)) , cx, cy) 



#9: 





1 


2 




f _pt ( ax 

i ay 

i bx i by , 

cx, cy) := 

(SOLUTI0NS(y = 

G(x , bx 

, by, cx , cy) a y = H(x, bx , by, cx, cy) 

, Lx, y 

D) 

#10: 









l 

Solv 

'ingy=G 

and 

y= H for (: 

x.y) gives 

the co-ordinate 

s of the 

Fermat point. 



Nov\ 

Net's try i 

with 

our spec! 

al values: 
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#11: f_pt(0, 0, 7, 2, 3, 5) 
#12: f_pt(0, 0, 7, 2, 3, 5) 


' 5731 2675-^3 3229-^3 

. 1226 3678 ' 3678 

[3.41483326, 3,129903399] 


1973 1 


1226 J 


We repeat the calculation for the Nspire-Triangle from page 19: 


#14: [2.344703633, 2,330623696] 

#15: triang1e(0, 0, 7, 2, 3, 5) 



The Cabri-plot with the rotated triangle ATB (which is used for the proof) 

Die Darstellungen zeigen deutlich, daB das Problem sicher nicht abschlieBend behandelt wurde. Es 
gibt noch viele Moglichkeiten zur Verfeinerung, zB beim Problem der Ableitungen oder bei der 
Schnittpunktsberechnung. Ich hoffe aber, daB deutlich wurde, wie ergiebig dieses Fermat-Problem fur 
den Mathematikunterricht in vielen Klassenstufen ist und wie hervorragend sich dabei ein Programm 
wie DERIVE einsetzen laBt. Der experimentelle Charakter der mathematischen Arbeit in unserem 
Projekt bewirkte jedenfalls einen erkennbaren Motivationsschub, der dem Mathematikunterricht ins- 
gesamt bestimmt zu gute kommt. 

(The presentations are showing very clear, that we couldn't deal with the problem suffi- 
ciently. There are many possibilities for refinement. But I hope that I could show how produc- 
tive this problem could be for many forms and how excellent you can use DERIVE for this 
purpose. The experimental character of the mathematical work in our project caused a sig- 
nificant motivation which could be very useful for maths teaching in total.) 


Additional comments for the revised version of DNL#16 will follow on the next page: 
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The students had the idea of applying calculus. So one could provide the suggestion to equate both 
partial derivatives to zero and . . . then try to solve the resulting system of equations: 


Inspecting the partial derivaties might discourage the students from solving the system. But they might 
also have heard or even learned to apply the Newton-Raphson-algorithm for approximatively solving 
equations. Why not extend for more variables?? 

Using the DERIVE-Help they might find NEWTONS (u , x , xO) . 


#32: Preci'si'onDi'gi' ts := 6 




Td d 



REVERSE 

NEWTONS 

— L(x, y), — L(x, y) 

1 [x. y]i [3, 3] 




..dx dy 

J 

JJ 


[3, 2, 1] 


#34: 


3 .41483 3.12990 " 

3.41483 3.12990 

3 .41483 3.12990 . 


Here we see the last three iterations of the algorithm. 

And finally it might be a welcome surprise to show the solution in a 3D-environment. 

L(x,y) from expression #6 is a function of two variables which can be plotted immediately in the 
3D-Plot window. The students have to find appropriate settings for the box. 

Using the Trace tool they are asked to find the “Minimum Point” and then they can add the point 

Ofp, Tfp? L(x F p, Tfp))- 



10.4 


z 
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Ebene Algebraische und 
Transzendente Kurven (5) 

Thomas Weth, Wurzburg, Germany 

Die Hundekurve - Die Konchoide des Nikomedes 

Die in dieser Folge behandelten Kurven sind Erweiterungen 
der sog. „Hundekurve“ oder „Tractrices“. Tractrices entste- 
hen, wenn ein Herrchen seinen Hund auf einem geraden Weg 
fiihrt, wahrend der Hund - bei konstanter Leinenlange - sei- 
nem „Lieblingsbaum“ zustrebt. Mathematisch formuliert 
handelt es sich bei dem Weg des Hundes um eine spezielle 
Form von Konchoiden. 

The curves dealt with in in this sequel are extension of the 
so called “Dog's curve“ or “Tractrices“. They are created 
when a man walks along a line while his dog - kept on a lead 
with constant length - wants to visit his favourite tree. From 
the mathematical point of view the dog's way is a special 
form of a conchoid.) 

Erzeugung von Konchoiden - How to create a conchoid 



Konchoiden im allgemeinen ergeben sich nach folgender 
Konstruktionsvorschrift: 

Gegeben ist eine ebene Kurve C und ein Punkt O in ihrer 
Ebene. Von einem Punkt Q der Kurve aus tragt man auf der 
Geraden OQ in die beiden moglichen Richtungen jeweils 
eine Strecke konstanter Lange K ab; die Endpunkte Pi und 
P 2 dieser Strecken sind dann Konchoidenpunkte zur Kurve 
C. 


This is the instruction how to obtain a conchoid: given is a plane curve and a point O in its plane. Have 
any point QeC. The two points P 1 and P 2 on the line OQ with QP 1 = QP 2 ’ = K (const) are points of a 
conchoid to C. If the curve C is a line the conchoids are especially simple, they are tractrices. Niko- 
medes (200 BC) invented this curve for doubling a die and for trisecting an angle. Its name is derived 
from the greek word for shell: ko jxt[.) 

Die Konchoide des Nikomedes - Herleitung der algebraischen Kurvengleichung 

Besonders einfache Konchoiden ergeben sich, wenn die Leitkurve C eine Gerade ist (wie im Beispiel 
der Tractrix). Diese Kurven wurden von Nikomedes (ca 200 v.Chr.) zur Wiirfelverdoppelung und zur 
Dreiteilung des Winkels erfunden. Der Name ist wegen der Kurvenform abgeleitet vom griechischen 
Koyxp, Muschel. 

Fur die Kurvenpunkte P 1/2 ergibt sich in Polarform (vgl. obige Zeichnung): 
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P • r = - 

^ 1,2 f 


-±K oder 


cos (p 


a 


\ 


f 


\ cos cp ) 


+ K = 0 und daraus fur die gesamte Kurve: 


r — - 


— K 

\ COS (p J 


\ 


r — - 


- + K 
y cos cp j 


= 0 oder kurz 


r — - 


= K 2 . 


y cos cp j 


(Polargleichung der Konchoide des Nikomedes; this is the polar form of Nikomedes' con- 
choid. Now you easily can find the Cartesian form:) 

Den Ubergang zu kartesischen Koordinaten erhalt man mit x-r cos (cp) und y-r sin(<^: 


f a-r V 
r 

V x ) 


= K 2 . Ersetzt man nun noch r 2 = x 2 + y 2 (DERIVE!) und faktorisiert die Summe, so 


erhalt man die algebraische Kurvengleichung 


{^x 2 + y 2 ^x-a ) 2 = K 2 x 2 . 

Aus der Kurvengleichung lafit sich nun entnehmen, dass die Kurven (fur verschiedene K) 

• symmetrisch zur x-Achse sind, 

• die Leitgerade C mit der Gleichung x = a zur Asymptote haben und 

• im Ursprung einen (evtl isolierter) Kurvenpunkten haben. 

(From the curve's equation you can see that for different values for K, the curves 

• are symmetric in respect to the x-axis, 

• have line C with x -a as an asymptote and 

• the origin is a (eventually isolated) point of the curve.) 


Fiir die Falle K< a, K = a und K> a ergeben sich die drei Kurvenformen: 



K< a K= a K> a 


DERIVE 6 has two remarkable features which we can use to reproduce the construction of 
this curve: the slider bars and the possibility to switch between rectangle coordinates and 
polar coordinates in one session. This one more occasion to narrow the gap between com- 
puter algebra and dynamic geometry. I'll try to explain the procedure very briefly guided by 
the following screen shots, Josef 
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Start in Rectangular coordinates. 

Plot #1 and #2 and introduce slider bars for a (0 < a < 10 ) 
and p_ (-tt/ 2 < p_< ir/2). 

#1: x = a 

#2: x-TANCpJ 

Moving the slider bars you can change a and see the rotating line. 
Change now to polar coordinates, 
introduce a third slider bar for k (0 < k < 1 0). 

Plotting #3 you can move the points on the two branches of the curve, 
exchanging p_by $ gives the polar form of the curve, 
and so does the implicit form in #4. 


You can alter a, k and p_ as you like. 

+ k p_ 


#3: 


C0S(p_) 


C0S(p_) 


k P- 


#4: r 


a V 2 
= k 


cosO) 



This is the “Dog's curve” on the hand held Tl-Nspire. 


( 1 .1 | 1 .2 | 1 .3 | DEG AUTO REAL □ 

7.5 

y\ 

\ k 
1 x 

- 10.36 \ 

k \ 

1 1 12.21 



♦ BSk 77 M= [i" 


( 1 .1 | 1 .2 | 1 .3 | DEG AUTO REAL «n H 

7.5 

y \ 

1 x 

k \ 

i 12.21 




^ 1:5, z 
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Verwendung der Konchoide zur Wiirfelverdoppelung - Doubling a die 

Zu konstruieren ist die Kantenlange x eines Wiirfels, der den doppelten Inhalt besitzt wie ein Wiirfel 
mit der Kantenlange a ; also x 3 = 2 a oder x = a 3 V2. 


Dazu zeichnet man eine Konchoide mit K = a und 
zeichnet durch U unter 60° eine Gerade zur x- 
Achse, die den rechten Kurvenzweig in S schnei- 
det. Die Gerade OS schneidet die Leitkurve C dann 
im Punkt T. x = OT ist dann die gesuchte Wiirfel- 
kante mit x 3 = = 2a. 

OT = x is the edge of the cube with x 3 = 2 a 3 . 

Nach dem Strahlenlsatz ist namlich 

x\a-2a\y bzw x-y = 2a 2 . (I) 

Mit dem Lot von O auf die Gerade US ist das 
Dreieck OSL rechtwinkelig und es gilt zunachst: 

LU = a cos (60°) = | . 



Damit folgt fur das Dreieck AOSL: (x + af = + OL 2 

und fur das Dreieck AOLU: a 2 = + OL 2 . 

4 

Subtraktion liefert: ( x + af - a 2 = + ~ f bzw. x 2 + lax = y 2 + ay 

(I) und (II) vereinfacht man durch Elimination von y zu x 3 = 2 a 3 . 


You can simplify (I) and (II) by eliminating y and you obtain x 3 = 2a 3 


Verwendung der Konchoide zur Drittelung des Winkels 

How to use the conchoid to trisect an angle 

1st a der zu drittelnde Winkel, so tragt man auf dem einen Schenkel eine beliebige Strecke [OU] mit 
der Lange a ab. Das Lot in U auf die Gerade OU schneidet den anderen Schenkel im Punkt N. 

Mit ON = c zeichnet man zur Leitgeraden NU eine Konchoide mit O als Pol und K = 2c. Durch N 

I I OL 

zeichnet man die Parallele zu OU, welche die Konchoide im Punkt P trifft. Dann ist: ZPOU = — . 
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If you want to trisect a, then take any arbitrary point U with OU = a on one side. The perpen- 
dicular intersects the other side of a in N. With ON = c draw a conchoid with respect to the 
line NU with O as ist pole and K = 2c. Let NP || OU , then ZPOU is the third of ZNOU 


Begrundung: 

Ist M die Mitte von QP, so liegen P, Q und N auf dem Thaleskreis iiber der Strecke [PQ]. Also ist 
PM = MQ = MN = c = ON. Damit gilt: 

I Z NPM |=| Z MOU I = fi (als „Z - Winkel“). 

Da das Dreieck MPN gleichschenkelig ist, gilt fur den AuBenwinkel ZOMN: 

| ZOMN| = 2 fi . 

Da auBerdem das Dreieck OMN gleichschenkelig ist, gilt wegen der Kongruenz der Basiswinkel: 


(M is in the middle of PQ, then P, Q and N are lying on a Thales circle, so PM = MQ = MN = 
c = ON. 

Then you will find two isosceles triangles AMPN and OMN A, giving simple relations between 
the angles.) 

| ZMON| = I Z OMN | = 2 fi. also ist a = 2 fi+fi oder fi = a/3. 


Es sei hier ausdriicklich darauf hingewiesen, daB es sich bei den Konstruktionen zur Wiirfelverdoppe- 
lung und zur Winkeldreiteilung nicht um elementare Zirkel- und Linealkonstruktionen im Pla- 
ton'schen Sinne handelt. Wie man an den vorhergehenden Beispielen sieht, ist namlich die Wiirfelver- 
doppelung und die Winkeldreiteilung sehr wohl geometrisch moglich - nur eben nicht - wie zB in der 
Galois-Theorie bewiesen wird - mit Zirkel und Lineal alleine. In den obigen Konstruktionen wurde 
namlich eine Gerade mit einer Kurve (der Konchoide) geschnitten, also eine Operation durchgefuhrt, 
die mit Zirkel und Lineal alleine nicht moglich ist. 


AbschlieBend sei bemerkt, dass sich fur einen Kreis als Leitlinie andere bekannte Konchoiden, u.a. die 
Pascalschen Schnecken und die Kardioide ergeben. Wir werden in einer spateren Folge darauf zu 
sprechen kommen. 


I want to emphasize that these constructions are not elementary constructions using ruler 
and compasses only in Platon s sense, because you need the intersection point between a 
conchoide and a line and this operation can 't be done with the tools mentioned. The Galois 
Theory proofs this fact. 


Finally I would like to mention, that using a circle as leading line we will find other well known 
conchoids as Pascal's snail and the cardiode. In a later sequel we will talk about. 
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DERIVE and plotting 
T - periodic functions 

Drs. J.C.M. Verhoosel, Eindhoven, Netherland 

Introduction. 


For most simple periodic functions it is rather easy to plot them. Think of y = sin(x) for instance. Just 
Author sin(x), press twice on Plot and voila. 

It is more difficult to plot functions like those in figure 1. 



We will discuss the possibilities and problems connected with plotting this kind of functions. 

The tools. 

Derive knows the following functions: 

a) CHI (a, x, b) ; a function which returns 1 for x e [a,b] and 0 outside [a, b] 

b) The IF - function. We can use this function for piecewise definded expressions using conditions. 
For instance: I F ( x < 0, x 2 ', x)isa function which for x < 0 returns a parabola and for x > 0 
the straight line y = x. 

c) MOD ( a , b ) which approximates to a modulo b. When MOD ( x , 1 ) is plotted, you get a saw- 
shaped- function. MOD (x, n) also returns a saw shaped function but now with period n. 

d) FLOOR (x) is narrowly connected to MOD (x) . Simplifying MOD (x, 1 ) yields x - FLOOR (x) . 
In other words: MOD (x, 1) = x - FLOOR (x). 

Plotting using CHI, IF and MOD 

The usage of CHI. 

Let's look at the following piecewise defined function: 

/(x) =x 2 ifx g [0,1] and 0 ifx g [0,1] 

We could Author this with: f(x) := CHI ( 0 , x , 1 ) x 2 . 

A short explanation: 

The fact that the product is forced to 0 outside [0,1] makes it rather difficult for us to plot functions 
only defined on an interval: DERIVE plots 0 outside the interval while the function is not defined over 
there. 
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Let's suppose we want a function like the one in figure 1. We can define J{x) on [0,1], on [1,2] etc. The 
next step is to define this function in terms of SUM: Apparently it's not a problem when the number 
of terms is finite. 


On [0,1] we get 4 CHI (0, x, 1) (x - 0.5) A 2. 

The quadratic function defined on [1,2] would be f (x) := 4 CHI(l,x,2) (x - 1.5) A 2; 

the one on [£, £+1]: 4 CHI(k,x,k+l) (x - 0.5-k) A 2. 

So when we take the first 5 terms we get: 

/(x) : =^ 4 CHI(£,x,£ + 1) (x - 0.5 - k) 2 

k=0 

This function isn't a true periodic function. To produce a true periodic function /(*) must become 

symmetric and extended to infinity (and negative infinity). So let us define g(x): 

00 

g(x) := ^ 4 CHl(k,x,k + 1) (x - 0.5 - k) 2 

k =—oo 

Although DERIVE is capable of plotting /(x), it accepts g(x) but seems not capable of plotting it. 
Probably DERIVE tries to evaluate the function and can't handle infinity. 


Using IF 

What about the IF - function? We can use this function to make a piecewise defined function: for in- 
stance: f (x) := IF(x < 0 OR x > l,x,x 2 ). 

A short explanation: 

When x < 0 or x > 1 , then f(x) = x, everywhere else/(*) = * 2 - 
Another example: g (x) := IF(x < -1 AND x < 2,x,?) 

When we plot g(x) we can see that outside [-1,2] g(x) is not forced to 0 because the function is not 
defined there. 

Nevertheless we have the same problems with the “IF - command 44 as those we encountered with the 
CHI-functions: it is not possible to plot true piecewise defined functions although it is possible to de- 
fine them! 


Using MOD(a^) 

Because MOD (x, a) already is a periodic function, it offers more perspectives. Moreover, it is not 
necessary for a to be an integer: for instance MOD (5,2.1) will yield 0.8. 

For students this function often gives surprises. It's not automatically clear for them why the plot of 
x MOD (x) or MOD (x A 2 ) looks so strange. These plots invite the students to discuss simplifications 
in DERIVE involving MOD (x) . Can we understand what is happening here? What is the connection 
between FLOOR (x) and MOD (x) ? After some experimenting in DERIVE they find it rather easy to 
understand that MOD (x) = x - FLOOR (x) and x MOD(x) simplifies to x A 2 - x FLOOR (x). 
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MOD(x) 


| 2D-plot 1:2 


x-FLOOR(x) 


plot 1 :4 


EB 

2D-plot 1:3 
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2D-plot 1:5 
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2D plot 1:6 


□00 


#1: 

MOD(x) 

\ y 

/ 

y 


#2: 
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;/ * 

WWWVv 

S/llll h 





#S: 

FLOOR (x) 



M0DO a 2) 


#6: 







Figure 2 


The usage of MOD(x) in plotting periodic functions. 

To understand how MOD(x,y) can be used to plot T-periodic functions, we must understand how 
MOD(x,y) works (see window 4 of figure 2 for a plot of MOD(x)). 

MOD(x) can be used to make sure that for instance/(3.2) will be reduced to/(0-2). This means that f(x) 
can be 1 -periodic: the part on the domain [0,1] is repeated. 

MOD(x) can be used to draw other, random composed, periodic functions. MOD(x) has period 1: 
every x is reduced to a value between 0 and 1. Replacing x with MOD(x) gives us a 1 -periodic func- 
tion! 

A nice way to create a T-periodic functions is to transform the function defined on [0,1] by multiplica- 
tion / translation etc. 

An example: 

We want a “periodic parabola 44 through the points (0,1) and (1,1) which is symmetrical on the axes 
x = 0.5. (see figure 3). 



Figure 3: the quadratic curve. 
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First we must find the formula of the quadratic curve: fix) := 4 (x - 0.5) 2 . This function surely is de- 
fined on [0,1]. Now we replace x by MOD(x) which gives us the function 
fix) := 4 (MOD(x) - 0.5) 2 . Plot the formula and you get figure 4. 



Suppose we want to translate this function so that the minimum of the quadratic curve occurs for 
x — 0, then we must replace x by x - 0.5. In the same way we can get a 2-periodic function by replac- 
ing x by x/2 etc. 


Another example: 

We want a series of half-circles. At first you can define a half circle on [0,1] (see figure 5). Solving for 
y gives us 2 possibilities, we choose the positive one. 

The only thing we must do now is to replace x by MOD(x). In the same way we get a 2-periodic series 
by replacing x by x/2. The period is twice as big but not the amplitude. So weTl have to multiply our 
function as well by a factor 2. Now we have a true 2-periodic series of half-circles. 



# 1 : 

# 2 : 
# 3 : 
#4: 
# 5 : 

# 0 : 

#7: 



2D-plot 1:1 


□0® 


2 


y 


3.5 ■ 


(x - 0.5) + y = 0.5 


2 2 2 

5QLVE((x - 0.5) + y = 0.5 , y) 


3 

2.5 


y = - i/Cx-(l - X)) V y = A*-C l - x)) 
Ax - (1 - x)) 

AMQD(x).(l - M0D(x))) 



.5 -3 - 2.5 -2 - 1.5 -1 - 0.5 


0.5 1 1.5 

- 0.5 


2.5 


3 . 



Figure 5: a series of half-circles 


(It is interesting that DERIVE “simplifies” the functions, when we have them annotated in the 
plot window! Josef) 


Noting that MOD(x,2) already has a period 2, it must be possible to do things more easy. 
Suppose we plot y = Vx Vl-x. 
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Suppose we plot y = 4x Vl-x. 

To make sure that the next period is produced we must be sure that x = 2. 1 for instance is reduced to 
x = 0. 1 . This can be done by MOD(x,2). So we replace x by MOD(x,2) and get: 


y = 7 MOD (x, 2) yjl - MOD (x, 2) . 


When we plot this function we get 2-periodic function in the next figure. 


Q 



1.5 




Let's try this on another function. Suppose/^) is defined as: 

/(x) = x 2 when x e [-0.2, 0.4] and 0.6 - periodic. 

We are tempted to replace x by MOD(x,0.6). Why does this not work? 


2 

M0D(x , 0.6) 



Looking at x = 0.7 we see that MOD(x,0.6) reduces it to 0.1 while it should be reduced to -0.1. We 
have to transform MOD(x). The period is all right: 0.6. The starting point is not correct. It should 
be -0.2, so: will MOD(x + 0.2,0. 6) work? Not yet. Suppose we look again at x = 0.7. Now it is re- 
duced to MOD(0.7 + 0.3, 0.6) = 0.3, still not correct. We must realize that we should subtract 0.2, so it 
should be: MOD(x + 0.2, 0.6) - 0.2! ! ! 

We can generalize this to the procedure shown in the next figure: 

2 

(M0D(x + 0.2, 0.6) - 0.2) 



Combining .... 

There is no restriction on the form of the function 
which is to be made periodic. We can join x 2 on 
[-0.2, 0.4] with 0.56 -x on [0.4,0. 6] for instance. 
In this case we make use of the IF- or CHI- 
function: 



■ 0.2 


■ 0.15 /\ 


■ 0.1 ■ / • \ • 


■ 0.05 ■ \ 

- 0 . 2 - 0.1 

0.1 0.2 0.3 0.4 0.5 \.6 


■- 0.05 
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We can generalize this to the procedure shown in the next figure: 

2 

x -x(-0.2, x, 0.4) + x(0.4, x, 0.6)-(-x + 0.56) 

Substitute x by MOD(x+0.2,0.8)-0.2 and plot the resulting expression! 

2 

(M0D(x + 0.2, 0.8) - 0.2) »xC-0 . 2 , M0D(x + 0.2, 0.8) - 0.2, 0.4) + x(0.4, 
MOD(x + 0.2, 0.8) - 0.2, 0.6)-(- (MOD(x + 0.2, 0.8) - 0.2) + 0.56) 



This works also with an IF-construction instead of the CHI-function for the piecewise defined 
function. Josef 


2 

IF(-0.2 < x < 0.4, x , IF(x < 0.6, 0.56 - x, ?)) 

2 

IF(-0.2 < M0D(x + 0.2, 0.8) - 0.2 < 0.4, (M0D(x + 0.2, 0.8) - 0.2) , 


IF(M0D(x + 0.2, 0.8) - 0.2 < 0.6, 0.56 - (M0D(x + 0.2, 0.8) - 0.2), ?)) 



How to automate the process. 

In lectures on Fourier analysis students construct a function /(0 by summing harmonics. 
When the original function is T-periodic, they can match it with their own constructed func- 
tion. To avoid unnecessary arithmetic one could automate the process by defining the next 
function: 

PERIODIC(f, t, a, b) := lim f 

t^M0D(t - a, b - a) + a 

From now on it's rather easy to use. 

PERIODIC(SIN(t), t, -1, 2) 

+ z 
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2 

PERIODIC (t , t, -0.5, 1) 



PERIODIC(EXPC-t), t, 0.3, tt) 



We have now the possibility to compare a Fourier - approximation of ’f[x) with/tr) itself, as 
you can see below! 

2 

F(t) := xC-1. t, l)-t + x(l. t, 2 . 5) • (t - 2 ) 

PERIODIC (F(t) , t, -1, 2.5) 



Conclusion. 

DERIVE can be a useful tool for plotting T-periodic functions. Although it is possible to define and 
use T-periodic functions with CHI and SUM they can't be plotted. The use of MOD(x,T) is simpler 
and more elegant. Nevertheless, the process of defining true T-periodic functions with CHI and SUM 
and the use of transformations is very leamingful. The use of PERIODIC! ) appears to be very useful 
in the classroom. 
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DERIVE, Newton Method and 
III - Conditioned Problems 

Jose M.Cardia Lopes, Porto, Portugal 

J. Bohm (1993) has recently written in this magazine about the Newton method for 
solving non linear equations pointing out, by means of an interesting example, the 
crucial importance of the choice of the initial values. The popularity of this method is 
due mostly to its quadratic convergence if the initial values are well chosen. (If it 
doesn't converge the divergence soon becomes obvious.) Nevertheless with ill- 
conditioned problems the Newton method (as any method ...) can lead to a large 
number of iterations and to not very good solutions. 

A particular class off ill-conditioned problems which often appear in the classroom is the prob- 
lem of solving polynomial equations with multiple roots. In fact it can be demonstrated that if we have 
an equation/f.v) = 0 with a multiple root X, if f p {x) is the computer representation of fix) (with the 
“perturbations” due to the ill-condition), if A is an upper bound on the perturbations in /(v), and if X’ is 
the corresponding root of f p {x) = 0, then 


\Xp -X\ 




where m is the multiplicity of the root. With this expression we can evaluate the effect of a perturba- 
tion in f{x) on the determination of the multiple root of the equation. As a result of the power 1/m the 
determination of multiple roots may be an ill-conditioned problem. 

Since the ill-conditioning is a consequence of the problem structure, we can only avoid it if we 
change the problem structure. For a polynomial equation we can substitute the evaluation of the multi- 
ple root of 

fix) = 0 

by the evaluation of the single root of 


where 


(see for example, Chapra e Canale, 1990). 


u(x) = 0 


h m = M 
{) /'(*) 


This is another situation in which DERIVE can be a very good pedagogical tool. We can 
graphically show how the structural change of the problem improves dramatically the path to the root. 
We can see some examples: 


Example 1 : /(x) = (x - 3) 6 (2x - 5) 


Derive has no difficulty to solve this equation but the graphical representation suggests that the New- 
ton method converges on the multiple root x = 3 at a very slow rate (notice that the function does not 
cross the x-axis at x = 3). The situation is very different with u(x). 
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Example 2 : f(x) = (5x - 1 4) 2 (x 2 + 3) 


Example 2: 

2 2 

#3: f(x) := (5-x - 14) ■ (x +3) 

f(x) 

u(x) := 

#4: d 

- f(x) 
dx 



If the lecture is designed for beginners there will probably be some interest in a more detailed 
approach (like the one prposed by Dyer (1993)) to show the results of some interations. 

In addition the second example points out the importance of the manipulation of the mathemati- 
cal expressions in order to minimize the number of arithmetic operations. See the plots of the factor- 
ized and the expanded function. 

(There is no difference now in higher versions of DERIVE, Josef) 



(Scale: x=5 10 -4 , y=5 10' 5 ; Accuracy: 9 ; Precision: Exact) 


Screen Shot from the 1994 DERIVE for DOS version 
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Exercise 


Find here a worksheet of Keith Eames: 


Use the graphical capability of Derive to answer these questions. Check your answers age- 
braically using Derive. 

Explain how you can get form F(x), by a single transformation whenever possible, to each 
G(x) in all of the following questions: 

1. Draw the graph of F(x) = x 2 - 2x + 3. 

a) G(x) = x 2 + 2x + 3 b) G(x) = - x 2 + 2x - 3 c) G(x) = - x 2 - 2x - 3 

d) G(x) = x 2 + 4x + 6 e) G(x) = 3x 2 - 6x + 9 c) G(x) = 4x 2 - 4x + 3 


2. Draw the graph of F(x) = (2x - 1) (3x + 5) 

a) G(x) = 6x2 + 7x + 4 b) G(x) = 6x 2 - lx - 5 c) G(x) = 6x 2 - 29x - 28 

d) G(x) = 10 - 14x - 12x 2 e) G(x) = 54x 2 - 2 lx - 5 f) G(x) = 6x 2 + 31x + 28 

3. Draw the graph of F(x) = 2x 3 - 3x 2 - 9x + 10 

a) G(x) = 2x 3 - 3x 2 - 9x b) G(x) = 2x 3 - 9x 2 + 3x + 14 

c) G(x) = - 4x J + 6x 2 + 18x - 20 d) G(x) = 10 + 18x - 12x 2 - 16x 3 

e) G(x) = 2x 3 + 9x 2 + 3x f) G(x) = 16x 3 - 12x 2 - 18x - 10 


4. Draw the graph of F(x) = sin (x) 

a) G(x) = 2 sin(x) cos(x) 

c) G(x) = cos 

e) G(x) = sin(x) + cos(x) 


b)G(x) = — sin(x) + -cos(x) 

d) G(x) = sin(x) - V2 cos(x) 

f) G(x) = 1+3 sin(x) - 4 sin 3 x 


All the trig functions in the above are in radians not degrees 


Hint: Use the well known fractional multiples of n 


5. Is it always possible to map a quadratic function to another quadratic function by a series of 
transformations? Experiment with your own quadratic equations. Can you prove any of 
your findings? Give some examples! What is the maximum number of different single 
transformations that is required? 
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Answers 

1. Use a trial and error method to find possible answers. 
Simplify F(x+1) 

Compare the given answers with G(x) 


a) F(-x) 

b) -F(x) 

c) - F(-x) 

2. If Derive does not give the answer shown 

a) F(x)+9 

b) F(-x) 

c) F(x-3) 

3. 

a) F(x) - 10 

b) F(x-l) 

c) - 2 F(x) 

4. 

a) F(2x) 

b) F(x + 71 / 6 ) 

c) F(2x) 


d) F(x+3) e) 3 F(x) 


d) - 2 F(x) e) F(- 3x) 


f) F(2x) 


f) F(x+2)-5 


d) F(-2x) 


e) F(x+2)+4 f) F(2x)-20 


d) 2 F(x - 7 i/ 4 ) e) V 2 F(x + tt/ 4) f) F(3x) + 1 


5. F(x) = ax 2 + bx + c can be transformed to G(x) = rx 2 + sx + 1 by at most 4 transformations, 
this depends on the values of a and r. 

If r and a are both positive then ist 3 transformations: 


s 2 b 2 

+ t — c — + 

4 r 4 a 



f 

b + s V 

v V a 

x — 

2 Vr Va 2 r J y 


How does this formula change if a is positive and r is negative? 

When you have finished this assignment and you are happy if your results and conclusions 
ask for the next worksheet. 


Answer the questions on it without using Derive. 


A correction: In the last DNL I published H. Scheuermann's phone number. One digit 
was wrong, this is the right number: 06192 27568 


Please don 't forget to renew your subscription for 1995 in time. Administration would be 
much easier for us. If you don ’ t cancel then we assume that you will stay a member. Many 
thanks for your cooperation. 

Bitte denken Sie daran, rechtzeitig Ihre Mitgliedschaft zu erneuern. Falls Sie nicht kiin- 
digen, diirfen wir annehmen, dass Sie Mitglied bleiben. Herzlichen Dank fur Ihr Ver- 
standnis. 


